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FEASIBILITT  OF  OPTIMIZING  NOZZLE  PERFORMANCE 
FOR  ORBITAL-LAUNCH  NUCLEAR  ROCKETS 
By  John  R.  Jack 


SUMMARY 

(  Nozzle  performance  for  orbital-launch  nuclear  rockets  was  evaluated  for  a 
range  of  operating  conditions  and  nozzle  geometries.  From  this  study,  conclu¬ 
sions  were  drawn  concerning  (l)  optimum  values  of  nozzle  area  ratio  and  operating 
pressure,  (2)  operating  conditions  for  favorable  nozzle  heat  transfer,  and 
(3)  the  effects  of  deviations  from  optimum  conditions.  ^ 

yjn  general,  the  nozzle  efficiency  increases  with  both  increasing  divergence 
angleand  increasing  stagnation-pressure  level.  For  a  given  set  of  conditions, 
however,  the  nozzle  efficiency  increases  with  increasing  area  ratio  only  until 
the  gains  in  performance  due  to  additional  expansion  are  offset  by  the  increase 
in  convective  heat  transfer.  The  variation  of  specific  impulse  with  pressure  is 
maximized  at  a  pressure  level  that  depends  upon  the  pressure  at  which  the  flow 
freezes.  On  the  basis  of  heat  transfer,  operation  on  the  low-pressure  side  of 
the  optimum  point  is  desirable.  This  mode  of  operation  sacrifices  very  little  in 
performance  but  permits  a  substantial  reduction  in  the  throat  heat  flux,  j 

INTRODUCTION 

y*A  basic  problem  associated  with  nuclear  rockets  is  the  optimization  of  the 
energy  absorbed  by  the  hydrogen  propellant  as  it  passes  through  the  reactor  and, 
in  turn,  the  maximization  of  the  power  available  for  thrust  (so-called  jet 
power )7j  Of  course,  the  amount  of  power  that  can  be  absorbed  by  the  hydrogen 
depends  upon  the  maximum  reactor  temperature  and  the  propellant  pressure.  Vith 
the  appropriate  reactor  operating  conditions  selected,  nozzle  operating  param¬ 
eters  are  desired  that  will  maximize  the  jet  power.  In  reference  l,Jit  was  dem¬ 
onstrated  that  moderate  changes  in  reactor  parameters  do  not  change  the  reactor 
performance  greatly  from  that  obtained  under  optimum  operating  conditions;  con¬ 
sequently,  nozzle-inlet  conditions  may  be  varied  somewhat  without  detracting  from 
the  reactor  performance.  Apparently,  sizable  dividends  may  be  available  by  op¬ 
timizing  nozzle  performance!"!  The  approach  to  be  considered  herein  will  be  to 
maximize  the  jet  power  and  then  to  determine  the  resulting  specific  impulse  and 
thrust. 

y At  the  high  temperature  levels  being  considered  (-<5000°  R)  for  nuclear  reac¬ 
tors,  two  major  nozzle  problems  are  encountered.  These  are  the  dissociation  of 
the  hydrogen  and  the  heat  transfer  to  the  nozzle  walls.  Both  of  these  problems 
are  intimately  related  to  the  temperature  and  the  pressure  of  the  propellant/j  IF 


the  flow  Is  assumed  to  be  in  thermodynamic  equilibrium,  a  reduction  in  operating 
pressure  at  constant  temperature  increases  the  enthalpy  level  and  consequently 
increases  the  specific  impulse.  Thus,  with  this  assumption,  it  is  advantageous 
to  operate  at  a  low  pressure,  which,  incidentally,  is  a  design  variable  subject 
to  the  designer's  choice.  On  the  other  hand,  if  the  flow  is  assumed  to  be  fro¬ 
zen,  it  is  desirable  to  operate  at  a  high  pressure  so  that  dissociation  is  inhib¬ 
ited  and  the  amount  of  energy  frozen  or  invested  in  dissociation  is  a  minimum. 
Actually,  however,  the  nozzle  flow  process  will  be  a  nonequilibrium  process,  so 
that  the  effect  of  pressure  on  nozzle  performance  is  not  well  defined  and  has  to 
be  considered  quite  carefully  in  any  analysis. 

The  need  for  considering  the  heat  transferred  to  the  nozzle  walls  in  a  noz¬ 
zle  study  is  well  exemplified  by  the  studies  of  reference  2.  The  heat-transfer 
analysis  of  reference  2  shows  that,  even  though  the  hydrogen  propellant  has  a 
large  heat  capacity,  a  nuclear-rocket  nozzle  cannot  be  regeneratively  cooled  for 
the  assumed  conditions  because  the  throat  heat  fluxes  encountered  are  so  high. 
Again,  the  indications  are  that  operation  at  a  low  pressure  level  is  desirable 
since  the  throat  heat  flux  would  be  low.  This  trend,  however,  cannot  be  fully 
accepted  until  all  important  nozzle  parameters  are  considered  together  and  worked 
into  a  nozzle  optimization  study. 

jjrhe  intent  of  this  report  is  to  make  a  general  nozzle  analysis  to  determine 
the  effect  of  frozen-flow  losses,  heaC^transfer  losses,  and  expansion  losses  on 
the  performance  characteristics  of  a  nuclear-rocket  nozzleT^n  particular^  nos?=^ 
zle  operating  conditions  and  geometry  will  be  so  varied  that  tentative  conclu¬ 
sions  concerning  optimum  operating  pressure  and  nozzle  geometry 
with  the  hope  that  better  performance  may  be  achieved  while  the 
transfer  is  reduced. 


ANALYSIS 

The  problem  to  be  considered  in  this  section  may  be  stated  as  follows:  The 
hydrogen  propellant,  upon  leaving  the  nuclear  reactor,  has  associated  with  it  a 
given  temperature  level  and  a  certain  amount  of  power.  What  nozzle  geometry  and 
operating  pressure  must  be  chosen  to  optimize  the  power  available  for  thrust  (jet 
power) ?  The  problem  may  be  defined  with  a  nozzle  power  balance. 


may  be  reached, 
convective  heat 


/ 


Power  Balance 


A  hypothetical  nuclear-rocket  nozzle  is  shown  in  the  following  sketch  with 
the  contributing  terms  of  the  power  balance  identified! 
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(All  symbols  are  conveniently  grouped  and  defined  in  the  appendix. )  The  propel¬ 
lant  (or  gas)  power  Pgas  must  he  divided  between  the  convective  power  loss 

Pq,  the  frozen  power  loss  P-p,  the  expansion  power  loss  PE,  and  the  jet  power 
Pv  Consequently,  the  nozzle  power  balance  is  given  by 

j 

Pgas  =  Pj  +  PF  +  PQ  +  PE  W 

An  overall  nozzle  efficiency  may  be  obtained  from  equation  (l)  by  forming  the 
ratio  of  jet  power  to  gas  power; 


PF  +  PQ  +  PE 
P 

gas 


(2) 


Equation  (2)  focuses  attention  upon  the  fact  that  the  various  nozzle  losses  must 
be  kept  to  a  minimum  in  order  to  achieve  maximum  nozzle  performance.  If  a  regen¬ 
erative  cooling  system  is  employed  to  minimize  the  convective  power  loss  Pq,  the 

reactor  power  required  to  yield  a  given  gas  power  is  less  by  the  amount  of 

power  recovered  regeneratively. 


Equation  (2)  may  be  cast  into  a  more  convenient  form  involving  a  frozen- flow 
efficiency  a  convective-heat-transfer  efficiency  t]q,  and  an  expansion  effi¬ 
ciency  Tjg : 


where  the  respective  efficiencies  are  given  by 

P  -  Rr, 

2 

gas 

..  (pgas  "  Pp)  “  PQ 

- P - I“b“ - 

gas  F 

(pgas  PF  Pq)  ~  PE  _ _ Pj 

%  '  Fgas  -  "  *gas  -  % 

The  problem  that  now  remains  is  to  find  P-g,  Pq,  and  Pg. 


-  P 


Q 


(3) 

(4) 

(5) 

(6) 


Frozen  Power  Loss 

During  the  expansion  of  the  hydrogen  (Hg)  through  the  rocket  nozzle,  atomic 
hydrogen  (H)  recombines  to  some  extent  and  thereby  releases  dissociation  energy 
and  increases  the  nozzle  performance  above  that  to  be  expected  without  recombina¬ 
tion  (i.e. ,  the  frozen-flow  case).  In  such  a  flow  situation,  the  performance 
gains  to  be  made  through  recombination  depend  to  a  great  extent  on  how  long  local 
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chemical  equilibrium  may  be  maintained#  Some  criterion  must  thus  be  employed  to 
predict  the  departure  of  the  nozzle  flow  from  chemical  equilibrium.  For  this 
purpose*  the  near-equilibrium  flow  criterion  proposed  in  reference  3  for  racket- 
nozzle  performance  calculations  will  be  used.  Specifically*  this  method  permits 
an  estimate  to  be  made  of  a  small  temperature  difference  representative  of  the 
difference  between  the  equilibrium  and  the  nonequilibrium  reactant  weight  frac¬ 
tions# 

The  reaction  of  interest  is 


% 

2H  +  M  ^  H2  +  M 


where  the  third  body  M  could  be  either  Hg  or  H.  The  relation  of  reference  3 
for  this  reaction  is  given  by 

^  2 


T*  -  T 


dT 

dt 


(is) 


i.  +  _i 
,%  % 


(7) 


where 


dissociation  rate,  liters^/(mole^) (sec) 
kr  recombination  rate,  liters2/(mole2) (sec) 

p  pressure,  atm 

R  universal  gas  constant,  (liters) (atm)/(°K) (mole) 

T  static  temperature,  °K 

T1  static  temperature  relating  actual  concentrations  to  equilibrium  values,  °K 
dT 

•Jk  rate  of  change  of  temperature  in  nozzle,  °K/sec 

Xj_  mole  fraction  of  constituent  i 


Rearranging  equation  (7)  and  assuming,  as  in  reference  3,  that  XH  <  0.1  and 
T*  _  T  20°  K  yield  the  following  equation  for  the  pressure i 


2.5 


(- 1)  (HT> 


20k. 


atm^ 


(8) 


Por  the  nuclear-rocket-nozzle  problem,  the  following  typical  parameters  have  been 
chosen  to  illustrate  the  approximate  pressure  level  required  to  assure  equilib¬ 
rium  flow: 


kr  1011  liters2/(mole2)(sec) 
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T  «  2800°  K 
-  «  10^°  K/sec 

at  1 

Substitution  of  these  parameters  into  equation  (8)  yields  p  >  l/3  atmosphere 
for  equilibrium  flow.  Actually  the  term  -dl/dt  varies  inversely  with  nozzle 
length;  considerable  latitude  in  the  pressure  required  to  assume  equilibrium  may 
thus  be  obtained  by  varying  the  nozzle  length.  The  use  of  the  criterion  of  ref¬ 
erence  3  does  not  imply  that  it  will  predict  the  point  of  departure  from  equilib¬ 
rium  accurately.  On  the  contrary,  at  best,  it  will  indicate  only  the  general 
region  of  the  flow  where  departures  from  equilibrium  are  likely.  For  the  pur¬ 
poses  of  this  report,  however,  the  method  should  be  adequate  to  demonstrate  the 
feasibility  of  optimizing  nozzle  performance. 

With  the  freezing  pressure  and  stagnation  conditions  known,  the  nozzle  pres¬ 
sure  ratio  at  which  freezing  occurs  can  be  found  and  with  the  aid  of  table  XI  of 
reference  4  the  degree  of  dissociation  a  can  be  calculated  from  the  following 
equation: 


a  - 


(9) 


Since  a  has  been  determined,  the  frozen  power  loss  may  be  determined  from 

PF  =  1.055  wHd,  kw  (10) 


where 

41,526a!!,- 

Hd  .  - jj - -<  Btu/lb 


(H) 


Convective -Heat- Transfer  Power  Loss 

A  nuclear-rocket  nozzle  is  heated  in  three  different  ways;  (l)  convective 
heat  transfer  from  the  propellant,  (2)  gamma  heating  from  the  reactor,  and 
(3)  thermal  radiation  from  the  reactor  face.  All  these  heat  losses  must  he  con¬ 
sidered  in  determining  the  total  nozzle  heat  load  (i. e. ,  the  total  amount  of  heat 
to  he  carried  away  hy  a  cooling  system)  hut  need  not  he  considered  in  defining 
nozzle  performance.  The  only  losses  that  must  he  considered  in  performance  cal¬ 
culations  are  those  that  stem  from  the  propellant,  that  is,  those  losses  that 
diminish  the  propellant  power  Pg^g.  Consequently,  the  heat  loss  of  interest 

herein  is  the  loss  associated  with  the  convective  heat  transfer  from  the  propel¬ 
lant  to  the  nozzle  walls. 

Ho  precise  calculating  procedure  (i. e. ,  one  that  accounts  for  dissociation* 
axial  pressure,  and  temperature  gradients)  exists  for  predicting  nozzle  heat 
transfer.  A  straightforward,  approximate  approach  that  has  been  used  success¬ 
fully  in  the  past  is  therefore  employed.  The  relation  used  to  predict  the  local 


5 


heat-transfer  coefficient  h  is  that  suggested  by  reference  5  for  fully  devel¬ 
oped  turbulent  pipe  flow: 


h  = 


0.023 

*0.2 


(0.067u  J°'2  .0.8 


ref' 

Pr2/3 

rrref 


/w\0,8/  T  \0*8 

\a)  (Trefj 


,  lb/(sq  ft) (sec) 


(12) 


where  Pref  and  Prref  are  evaluated  at  a  reference  temperature  Tref  defined 
as 


Tref  »  \  (*  +  V  (13) 

The  convective  heat  flux  may  then  be  obtained  from 

q  «  h(H0  -  H^),  Btu/(sq  ft) (sec)  (14) 


The  total  nozzle  heat  loss  that  is  required  for  evaluating  nozzle  perform¬ 
ance  may  be  estimated  conveniently  by  the  method  of  reference  6,  in  which  the 
total  convective  power  loss  to  a  conical  nozzle  in  terms  of  the  local  heat  flux 
at  the  nozzle  throat  is  given  by 


In" 


PQ  a  1.055  qtAt  0.8  +  — - 


lew 


(15) 


The  use  of  equation  (15)  simplifies  the  calculation  considerably  since  only  the 
throat  heat  flux  need  be  determined.  The  remaining  parameters  are  specified 
through  the  nozzle  geometry. 


Expansion  Power  Loss 

The  expansion  power  loss  is  determined  by  the  amount  of  thermal  energy  re¬ 
maining  in  the  propellant  flow  at  the  exit  plane  of  the  nozzle.  This  power  loss 
is  thus  directly  related  to  the  expansion  area  ratio/  the  greater  the  expansion 
area  ratio,  the  smaller  the  expansion  power  loss.  There  is  no  restriction  on  the 
expansion  area  ratio  due  to  ambient  pressure  as  this  pressure  is  taken  equal  to 
zero  in  orbit. 


Calculations  of  these  losses  assume  a  one-dimensional  isentropic  expansion 
of  the  propellant  through  the  nozzle.  In  reality,  the  isentropic  assumption  is 
in  conflict  with  the  previous  discussion  concerning  convective-heat-transfer 
losses.  Even  so,  reference  3  indicates  that  heat  transfer  to  the  nozzle  walls 
introduce  only  small  errors  into  the  calculated  performance  for  rocket  en¬ 
gines  of  reasonable  size.  The  expansion  power  loss  is  found  from  the  following 
equation: 

* 

P-g  sa  1.055  wH-g,  kw 

vhere  H*g  is  determined  by  nozzle  stagnation  conditions  and  nozzle  expansion 
area  ratio. 
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Optimization  Procedure 


The  performance  calculations  upon  -which  the  optimizations  and  comparisons 
are  based  were  made  by  specifying  the  following  basic  initial  parameters: 
reactor- outlet  temperature  TQ,  reactor  power  level  Pgas,  and  reactor-outlet 

pressure  level  pQ.  The  variables  to  be  optimized  are  expansion  area  ratio 
Ae/A^  nozzle  efficiency  r\  or  jet  power  Pj,  and  finally  the  stagnation  pres¬ 
sure  pQ. 

The  method  will  now  be  outlined  briefly*  With  the  reactor- outlet  tempera¬ 
ture  TQ  and  a  selected  stagnation  pressure  pQ>  the  stagnation- enthalpy  level 

H0  can  be  determined  from  reference  4*  Now,  with  the  selected  operating  condi¬ 
tions  and  the  criterion  of  reference  3,  the  approximate  pressure  at  which  the 
flow  freezes  can  be  found.  This  pressure  then  determines  the  nozzle  area  ratio 
(see  ref*  4)  at  which  the  flow  freezes  and  permits  a  to  be  determined 
(eq.  (9))*  Once  the  "freezing"  area  ratio  is‘  known,  the  static- enthalpy  distri¬ 
bution  through  the  nozzle  may  be  calculated.  With  this  information^  the  frozen 
power  loss  may  be  found  from  equation  (10) ,  and  the  expansion  power  loss  at  any 
desired  exit  area  ratio  may  be  determined  from  equation  (16)* 


All  that  remains  to  be  found  in  order  to  determine  nozzle  performance  is  the 
power  invested  in  convective  heat  transfer  Pq.  Finding  this  power  loss  requires 

the  nozzle  flow  rate  and  the  throat  area*  The  nozzle  flow  rate  is  given  by  the 
reactor  power  level  and  the  stagnation  enthalpy: 


W  a 


P 

gas 

1.055  H,/ 


lb/sec 


(17) 


If  the  flow  through  the  nozzle  throat  is  in  equilibrium,  as  determined  by  the 
criterion  of  reference  3,  the  throat  area  is  found  by  calculating  an  isentropic 
expansion  from  the  stagnation  conditions  and  maximizing  the  product  pu.  The 
throat  area  is  then  given  by 

A+  =  — ,  sq  ft  (18) 

0  pu 

For  the  case  in  which  the  flow  through  the  throat  is  frozen,  the  throat  area  is 
determined  from  the  perfect-gas  equation 


where  Z  =  1  +  a  has  been  inserted  to  account  for  compressibility  and  both  a 
and  r  are  determined  by  the  conditions  at  the  freezing  area  ratio.  With  this 
information,  the  convective  power  loss  Pq  is  readily  found  from  equations  (12) 

to  (15). 
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After  all  the  powers  required  to  determine  performance  have  been  found,  the 

specific  impulse  and  the  thrust  may  be  found  from  the  jet  power  p.  as  follows: 

J 

Pj  “  ^gas 


and 


TIP 


gas 


sec 


w 

P  fa  wl,  lb 


(20) 


Equations  (20)  neglect  the  pressure-area  term  in  the  momentum  equation;  this  is 
justified,  however,  when  it  is  recalled  that  the  ambient  pressure  is  taken  equal 
to  zero  in  a  space  orbit  and  that  the  nozzle-exit  Mach  number  is  large,  that  is, 
the  nozzle- exit  static  pressure  is  quite  small. 


By  use  of  the  procedure  just  outlined  and  the  assumptions  and  relations  de¬ 
scribed  in  the  previous  sections,  nozzle  efficiencies  Were  calculated  for  the 
desired  combination  of  independent  variables  and  the  maximum  nozzle  efficiencies 
were  found. 


The  following  values  and  ranges  of  variables  were  used: 


Reactor- outlet  temperature,  T  ,  °R 
Reactor- outlet  pressure,  pQ,  atm  .  . 
Reactor  power  level,  Pgas,  kw  .  .  , 
Conical  nozzle  half-angles,  e  ,  deg 
Nozzle  wall  temperature,  1^,  °R  .  , 


.  5040 

.  0.01  to  100 

*  *  *  •  •  7 000 
10,  15,  20,  25 

•  •  *  •  •  14:40 


/ 

These  parameters  are  typical  of  those  chosen  for  the  low-power  space  nuclear 
rocket  discussed  in  reference  7. 


RESULTS  AND  DISCUSSION 

_  ®ie  Primary  purpose  of  this  analysis  was  to  determine  the  nozzle  operating 
conditions  and  the  geometry  that  would  produce  the  optimum  nozzle  performance. 
Emphasis  was  also  placed  on  the  total  convective  heat  load  and  the  nozzle- throat 
heat  flux,  since  both  may  constitute  a  major  design  limitation. 

As  might  be  surmised  from  the  ANALYSIS,  the  effect  of  freezing  pressure  upon 
nozzle  performance  is  quite  important.  For  a  better  understanding  of  the  nozzle 
theref°re;>  the  discussion  is  initiated  with  consideration  of  the  general 
effect  of  freezing  pressure.  Then,  for  a  detailed  discussion  of  nozzle  perform- 
ance,  a  specific  but  realistic  freezing  pressure  (based  on  the  best  data  avail¬ 
able)  is  chosen  to  make  complete  nozzle  performance  calculations. 
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Effect  of  Freezing-Pressure  Level 

According  to  the  near- equilibrium  criterion  of  reference  3,  a  local  nozzle 
pressure  greater  than,  or  equal  to,  a  certain  value  (Pfr)  must  be  obtained  in 

order  for  equilibrium  flow  to  exist.  If  the  local  pressure  falls  below  pfr,  the 

flow  departs  from  equilibrium,  approaches  a  frozen  condition  in  a  short  transi¬ 
tional  region,  and  finally  is  completely  frozen.  For  analysis  purposes,  the  flow 
will  be  assumed  to  freeze  suddenly  at  p^  .  Thus,  if  the  nozzle-inlet  pressure 

is  greater  than  p_^,  th e  flow  will  be  in  equilibrium  in  the  nozzle  until  it  ex¬ 
pands  to  a  pressure  less  than  or  equal  to  *  The  expansion  from  this  point 
on  is  then  completely  frozen.  If  the  nozzle- inlet  pressure  is  less  than  p^  , 

the  flow  is  completely  frozen.  The  pressure  at  which  the  flow  freezes  should 
thus  affect  the  nozzle  performance  considerably.  As  the  inlet  pressure  increases 
and  exceeds  the  freezing  pressure,  the  overall  nozzle  efficiency  should  depart 
from  the  frozen  curve  and  approach  the  equilibrium  curve. 

The  effect  of  freezing  pressure  may  be  observed  in  figure  1  for  a  nozzle 
half-angle  0C  of  15°,  an  area  ratio  Ae/A^  of  50,  and  for  three  freezing  pres¬ 
sures:  <0.01,  l/3,  and  100  atmospheres.  For  some  freezing  pressure  less  than 
0.01  atmosphere,  the  nozzle  flow  will  be  in  complete  equilibrium  for  the  inlet 
pressures  considered.  (This  curve  is  presented  as  a  reference  curve  and  repre¬ 
sents  the  best  achievable  performance  for  the  conditions  investigated. )  For  a 
pfr  of  100  atmospheres,  the  flow  will  be  completely  frozen,  since  all  inlet 

pressures  considered  are  less  than  100  atmospheres  and  consequently  all  nozzle 
pressures  will  also  be  less  than  100  atmospheres.  Also  presented  in  figure  1  is 
the  nozzle- efficiency  variation  for  the  freezing  pressure  (l/3  atm)  derived  with 
the  aid  of  equation  (8).  As  the  pressure  increases  and  exceeds  l/3  atmosphere, 
the  flow  departs  from  the  completely  frozen  curve  and  the  nozzle  efficiency 
approaches  the  equilibrium  efficiency.  The  equilibrium  curve  is  reached  at  a 
pressure  level  of  approximately  3  atmospheres,  and  then  the  two  curves  are  the 
same.  As  will  be  shown  later,  operation  at  a  low  pressure  has  the  advantage  of 
decreasing  the  throat  heat  flux  greatly,  which  thereby  helps  to  alleviate  a 
heat-transfer  problem.  On  the  other  hand,  the  total  heat  load  increases j  how¬ 
ever,  this  increase  poses  no  real  design  problem.  All  the  following  results  to 
be  presented  are  based  on  the  value  of  p^r  of  l/3  atmosphere,  which  was  derived 
in  the  ANALYSIS. 


Effect  of  Expansion  Area  Ratio 

The  effect  of  nozzle  expansion  area  ratio  on  overall  nozzle  efficiency  for 
a  given  stagnation  pressure  and  divergence  angle  is  shown  in  figure  2  for  stag¬ 
nation  pressures  of  0.01,  1,  and  100  atmospheres.  In  general,  the  nozzle  effi¬ 
ciency  increases  as  the  expansion  area  ratio  increases.  A  point  is  finally 
reached,  however,  at  which  an  additional  increase  in  area  ratio  decreases  the 
nozzle  efficiency.  At  this  point,  the  gain  in  performance  from  additional  ex¬ 
pansion  is  offset  by  the  increase  in  nozzle  convective  heat  loss.  In  addition, 
figure  2  indicates  that  the  nozzle  efficiency  increases  with  increasing  nozzle 
divergence  angle  6Q  and  that  the  optimum  efficiency  at  a  given  6C  progresses 
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to  higher  area  ratios  as  the  divergence  angle  increases.  Even  though  maximum 
nozzle  efficiencies  have  been  found,  the  penalty  associated  with  moderate  devia¬ 
tions  from  the  maximum  is  not  large,  since  each  curve  is  rather  flat  in  the  re¬ 
gion  of  the  maximum  efficiency.  An  area  ratio  less  than  optimum  could  therefore 
be  chosen  to  reduce  nozzle  weight. 


Effect  of  Stagnation-Pressure  Level 

The  variation  of  maximum  nozzle  efficiency  from  figure  2  with  stagnation- 
pressure  level  is  illustrated  in  figure  3.  Eor  a  given  divergence  angle,  an  in¬ 
crease  in  stagnation  pressure  produces  a  substantial  increase  in  nozzle  effi¬ 
ciency.  A  stagnation-pressure  level  is  finally  reached,  however,  after  which  an 
increase  in  pressure  produces  a  very  small  increase  in  nozzle  efficiency.  Thus, 
for  the  ease  being  considered,  gains  in  performance  are  negligible  for  stagnation 
pressures  greater  than  about  10  times  the  freezing  pressure.  The  same  behavior 
pattern  occurs  for  all  divergence  angles,  the  only  difference  being  that  the 
absolute  level  for  the  nozzle-efficiency  curve  increases  as  the  divergence  angle 
increases.  At  a  pressure  of  3  atmospheres,  changing  the  divergence  angle  from 
10°  to  25°,  for  example,  changes  the  maximum  nozzle  efficiency  from  0.85  to  0.92. 

The  change  in  thrust  level  at  optimum  nozzle  efficiency  with  stagnation 
pressure  is  shown  in  figure  4.  The  preceding  discussion  concerning  nozzle  effi¬ 
ciency  is  directly  applicable.  The  only  additional  point  to  be  made  is  that  the 
optimum  thrust  level  does  not  change  much  with  nozzle  divergence  angle. 

The  most  interesting  effect  of  stagnation-pressure  level  may  be  found  in 
figure  5,  in  which  the  variation  of  specific  impulse  with  pressure  at  maximum 
nozzle  efficiency  is  presented  for  divergence  angles  of  10°  and  25°.  The  spe¬ 
cific  impulse  is  optimized  at  a  pressure  level  of  approximately  3  atmospheres  for 
both  divergence  angles.  This  result,  of  course,  is  a  consequence  of  the  freezing 
pressure  used  in  the  analysis.  The  optimum  pressure  level  would  vary  with  the 
freezing  pressure. 

The  previous  discussion  has  been  concerned  only  with  the  optimum  parameters. 
This  approach  is  not  too  realistic  because  the  optimum  area  ratio  increases 
greatly  as  the  stagnation  pressure  increases.  At  a  stagnation  pressure  of 
100  atmospheres,  for  example,  the  optimum  area  ratio  for  a  divergence  angle  of 
25°  is  of  the  order  of  20,000.  It  is  therefore  of  interest  to  consider  more 
practical  area  ratios,  as  in  figure  6,  in  which  specific  impulse  is  given  in 
terms  of  stagnation  pressure  and  area  ratio  for  a  divergence  angle  of  15°.  Con¬ 
sidering  area  ratios  other  than  the  optimum  does  not  change  the  character  of  the 
curves.  Reducing  the  area  ratio,  however,  slightly  reduces  the  maximum  specific 
impulse  obtainable.  At  a  stagnation  pressure  of  3  atmospheres,  for  example,  the 
optimum  specific  impulse  is  only  1.06  times  that  obtained  for  an  area  ratio  of 
25.  Apparently,  area  ratios  varying  all  the  way  from  the  optimum  value  to  25  may 
be  employed  without  a  significant  decrease  in  specific  impulse  (although  nozzle 
efficiency  decreases  considerably). 


Effect  of  Heat  Transfer 


As  noted  in  figure  5,  the  maximum  specific  impulse  is  not  much  greater  than 
those  obtained  at  pressures  greater  than  or  slightly  less  than  3  atmospheres. 
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’This  fact  is  significant  in  terms  of  the  nozzle  heat-transfer  problem  for  the 
following  reasons:  The  feasibility  of  nozzle  cooling  for  a  given  set  of  condi¬ 
tions  is  primarily  determined  by  the  throat  heat  flux,  since  it  is  the  highest 
heat  flux  experienced  by  a  nozzle*  If  regenerative  cooling  is  employed,  the 
coolant  flow  must  be  capable  of  handling  the  highest  nozzle  heat  flux,  namely, 
the  throat  heat  flux*  This  requirement  is  rather  severe,  as  noted  in  refer¬ 
ence  2,  in  which  this  particular  aspect  of  regenerative  cooling  was  investigated. 
If  radiation  cooling  is  employed,  the  local  heat  flux  must  be  radiated  away  at  a 
nozzle  wall  temperature  that  does  not  exceed  the  nozzle-wall-material  melting 
temperature.  These  considerations  are  the  most  important  ones/  a  secondary 
effect,  however,  that  must  be  considered  for  the  case  of  regenerative  cooling  is 
the  total  cooling  load*  As  noted  in  equation  (15),  this  effect  is  governed  also 
by  the  throat  heat  flux/  however,  since  hydrogen  has  a  good  heat  capacity,  this 
effect  is  not  so  important  as  the  local  heat  flux  at  the  nozzle  throat* 

As  the  stagnation  pressure  is  increased,  the  throat  heat  flux  increases/  in 
fact,  it  increases  approximately  as  the  0. 75  power  of  the  stagnation  pressure 
K  »  329  pg-75).  it  therefore  most  desirable  to  operate  at  the  lowest  pres¬ 
sure  possible  without  sacrificing  greatly  on  performance*  Figure  7  shows  that  by 
operating  at  the  pressure  level  yielding  the  maximum  specific  impulse  (3  atm), 
the  throat  heat  flux  may  be  reduced  by  over  a  factor  of  10  from  the  value  ob¬ 
tained  at  a  pressure  of  100  atmospheres.  In  addition,  by  operating  at  a  pressure 
slightly  lower  than  the  optimum  pressure,  the  specific  impulse  is  not  changed 
greatly,  but  the  throat  heat  flux  is  changed  considerably.  Operation  at  1  atmos¬ 
phere,  for  example,  reduces  the  specific  impulse  by  approximately  4  percent  and 
reduces  the  throat  heat  flux  by  a  factor  of  2  below  that  obtained  at  the  optimum 
stagnation  pressure.  A  small  sacrifice  in  specific  impulse  thus  reduces  the 
throat  heat  flux  by  a  large  amount.  Although  the  sacrifice  In  specific  impulse 
Is  not  large,  reference  to  figure  3  shows  that  the  optimum  efficiency  is  reduced 
by  approximately  13  percent,  which  requires  a  higher  reactor  power  for  a  fixed 
thrust  and  specific  Impulse. 

The  effect  of  stagnation-pressure  level  on  the  total  convective  cooling  load 
may  be  seen  In  figure  8  for  a  divergence  angle  of  15°  and  an  area  ratio  of  50. 
Reducing  the  pressure  increases  the  total  convective  heat  load  because  the  throat 
area  increases  faster  than  the  throat  heat  flux  decreases  (see  eq.  (15)).  This 
increase  in  Q,  however.  Is  slight  and  is  not  so  Important  a  consideration  as  the 
throat  heat  flux  q^. 

At  this  point  it  should  be  reemphasized  that  the  heat  load  presented  In  fig¬ 
ure  8  does  not  represent  the  total  heat  load  that  a  cooling  system  would  have  to 
handle.  The  heat  due  to  gamma  rays  and  thermal  radiation  from  the  face  of  the 
reactor  is  lacking.  Estimates  of  these  heat-transfer  effects  may  be  made,  how¬ 
ever.  In  reference  3,  it  is  shown  that  the  nuclear  heating  can  be  as  much  as 
10  percent  of  the  reactor  power  at  a  reactor  power  level  of  50  megawatts.  As 
the  reactor  power  decreases,  however,  the  nuclear  heating  also  decreases.  For 
the  present  case  (7000  kw),  it  is  therefore  estimated  that  the  nuclear  heating  is 
negligible  compared  with  the  other  heating  effects.  The  term  that  remains  to  be 
estimated  is  that  due  to  thermal  radiation.  This  term  has  been  estimated  by 
utilizing  the  data  of  reference  2,  and  its  value  is  approximately  440  Btu  per 
second.  Since  the  thermal  radiation  effect  is  not  a  function  of  the  pressure 
level,  440  Btu  per  second  may  be  added  to  the  curve  in  figure  8  to  obtain  an 
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estimate  of  the  total  heat  load  to  he  handled  hy  a  cooling  system  for  the  condi¬ 
tions  investigated  herein. 


Effect  of  Power  Level 

The  immediate  effect  of  changing  the  power  level  for  a  given  reactor-outlet 
temperature  TQ,  reactor- outlet  pressure  pQ,  and  Wall  temperature  T^  is  to 

vary  the  propellant  flow  rate  w.  This  change  in  flow  rate,  in  turn,  affects 
the  nozzle  performance  through  the  convective  heat  transfer  and  the  thrust. 

The  new  flow  rate  w2  is  related  to  the  flow  rate  used  in  the  analysis  w-^ 

Ly 


gas,l 


With  the  use  of  equation  (21)  and 


h,  2  =  kAfc,l  (22) 

the  local  throat  heat  flux  and  the  total  nozzle  heat  load  are  found  from  equa¬ 
tions  (14)  and  (15)  to  he 


(23) 

(24) 

(25) 


Investigation  of  the  effect  of  k  on  r)^2  for  typical  values  of  ru  1  and  for 
values  of  k  ranging  from  1  to  100  shows  that  rn  _  does  not  vary  significantly 
from  t)q^,  a"k  most  about  6  percent.  As  the  propellant  power  is  increased,  the 

specific  impulse  therefore  remains  almost  constant  (as  would  be  expected  with  a 
constant  stagnation  temperature  Tq)  and  the  thrust  increases  directly  as  the 

power  increases;  that  is,  P2  =  kFq.  Thus,  the  maximum  specific  impulses  found 
at  a  Pgas  of  7000  kilowatts  should  also  be  approximately  right  for  higher  power 
levels. 


The  throat  heat  flux  is  reduced  as  the  power  level  is  increased.  With 
k  =  100,  for  example,  the  throat  heat  flux  is  0.63  of  the  value  obtained  for  a 
power  level  of  7000  kilowatts. 
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Effect  of  Nozzle-Wall  Temperature 

Since  the  nozzle-wall  temperature  has  been  considered  an  independent  design 
variable,  its  effect  upon  the  overall  nozzle  efficiency  should  be  considered. 
This  effect  will  not  be  treated  in  detail  but  in  a  general  manner,  and  some  gen¬ 
eral  conclusions  concerning  it  will  be  made. 

As  the  wall  temperature  is  increased,  the  driving  potential  for  the  local 
heat  flux  H0  -  is  reduced,  and  thereby  the  local  heating  rate  is  reduced. 
The  net  result  is  that  the  heat-transfer  efficiency  is  increased  with  a 

corresponding  increase  in  overall  nozzle  efficiency.  The  total  amount  of  power 
transferred  to  a  hot  wall  in  terms  of  the  power  transferred  to  a  cold  wall  is 
given  approximately  by 


Q,hot 


(H0 


Hw) 


hot 


v 


cold 


Q, cold 


This  expression  is  not  exact,  of  course,  since  the  heat-transfer  coefficient  h 
would  change  slightly  with  wall  temperature. 


CONCLUSIONS 

The  feasibility  of  optimizing  nozzle  performance  for  an  orbital-launch  nu¬ 
clear  rocket  was  considered.  Power  losses  due  to  dissociation,  heat  transfer 
from  the  propellant  to  the  engine  walls,  and  incomplete  expansion  of  the  propel¬ 
lant  were  evaluated.  Operating  conditions  used  (reactor  power  level,  7000  kw  and 
reactor- outlet  temperature,  5040°  R)  are  those  under  consideration  for  low-power 
space  nuclear  rockets.  In  addition,  a  partly  frozen  nozzle  flow  was  studied. 

The  following  conclusions  based  on  the  analysis  presented  herein  may  be  made: 

1.  The  nozzle  efficiency  for  a  given  nozzle  divergence  angle  and  pressure 
increases  with  increasing  area  ratio  until  the  gains  in  performance  due  to  addi¬ 
tional  expansion  are  offset  by  the  increase  in  convective  heat  transfer.  At  this 
area  ratio,  the  nozzle  efficiency  starts  to  decrease.  In  general,  the  nozzle 
efficiency  increases  with  both  increasing  divergence  angle  and  stagnation- 
pressure  level. 

2.  For  a  given  divergence  angle,  the  maximum  nozzle  efficiency  increases 
with  increasing  stagnation-pressure  level.  A  stagnation  pressure  is  finally 
reached,  however,  at  which  the  variation  of  nozzle  efficiency  with  pressure  is 
small.  The  optimum  thrust  level  behaves  in  a  similar  manner  with  the  stagnation- 
pressure  level. 

3.  The  specific  impulse  Is  maximized  at  a  pressure  level  that  depends  upon 
the  pressure  at  which  the  flow  freezes.  For  the  case  treated  herein,  in  which 
the  flow  freezes  at  a  pressure  of  l/3  atmosphere,  the  best  stagnation-pressure 
level  is  about  10  times  the  freezing  pressure.  On  the  basis  of  heat  transfer, 
operation  on  the  low-pressure  side  of  the  optimum  point  is  desirable.  This  mode 
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of  operation  sacrifices  very  little  in  specific  impulse  Tout  causes  a  substantial 
reduction  in  the  throat  heat  flux. 

4.  The  throat  heat  flux  increases  substantially  with  increasing  stagnation 
pressure  and  varies  approximately  as  the  0.75  power  of  the  stagnation  pressure. 

5.  The  total  heat  load  to  be  handled  by  a  cooling  system  increases  with 
decreasing  pressure. 

Lewis  Research  Center 

National  Aeronautics  and  Space  Administration 
Cleveland,  Ohio,  November  30,  1962 
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APPENDIX  -  SYMBOLS 

A  flow  area,  sq  ft 

d  diameter ^  ft 

Ej_  dissociation  potential,  v 
P  thrust,  lb 

g  acceleration  due  to  gravity,  ft/sec^ 

H  enthalpy ,  Btu/lb 

I  effective  specific  impulse,  sec 

M  molecular  weight,  lb 

P  power,  kw 

Pr  Prandtl  number 

p  pressure,  atm 

Q  total  heat  flux,  Btu/sec 
q  local  heat  flux,  Btu/(sq  ft) (sec) 

R  hydrogen  gas  constant 

T  temperature,  °R 

t  time,  sec 

u  stream  velocity,  ft/sec 

w  propellant  flow  rate,  lb/sec 
a  degree  of  dissociation 

X  ratio  of  specific  heats  for  frozen  flow 

T)  nozzle  efficiency 

6C  conical  nozzle  half-angle  (divergence  angle),  deg 
ji  absolute  viscosity,  micropoises 
p  stream  density,  Ib/cu  ft 

Subscripts: 
d  dissociation 


E  expansion 
e  nozzle  exit 

F  frozen 

fr  freezing 

gas  propellant 

j  jet 

o  reactor  outlet  (nozzle  stagnation) 

Q  heat  transfer 

ref  reference  condition  for  heat  transfer 
t  nozzle  throat 

v  nozzle  wall 
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Expansion  area  ratio,  Ae/At 

(c)  Stagnation  pressure,  100  atmospheres. 

Figure  2.  -  Variation  of  nozzle  efficiency  with  nozzle  geometry. 
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Figure  7.  -  Nozzle -throat  heat  flux 
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